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ABSTRACT 

Using the derived gamma-ray burst -E p0 ak and fluences from the complete BATSE 5B Spectral 
Catalog, we study the ensemble characteristics of the -E pe ak-fluence relation for GRBs. This 
relation appears to be a physically meaningful and insightful fundamental discriminator between 
long and short bursts. We discuss the results of the lower limit test of the E pea k-E- 1BO relations 
in the -Epeai^huence plane for BATSE bursts with no observed redshift. Our results confirm the 
presence of two GRB classes as well as heavily suggesting two different GRB progenitor types. 

Subject headings: gamma rays: bursts — methods: data analysis 



1. Introduction 

Classification of Gamma-Ray Bursts (GRBs) is 
certainly a difficult task. Bursts are divided into 
long and short classes, based upon the bimodal 



durat ion histogram of bursts (|Kouveliotou et al 



Il993l ) observed by the Burst And Transient Source 
Experiment (BATSE), which was on board the 
Compton Gamma Ray Observatory. One param- 
eter, the spli t time of 2 seconds on the £90 du- 
rations plot ( Koshut et al. 19961 ). was sufficient 



to classify bursts. Bursts from the two classes 
have another discriminator, spectral hardness as 
determined by the ratio of two broad energy chan- 
nels. This hardness ratio, when used in conjunc- 
tion with the duration, provides a means for clas - 
sification, as shown by iKouveliotou et al I (|l993h . 
In addition, another classification scheme uses 
the scatter plot of the fluen ce and durat i on fit - 
ted with two 2D Gaussians (jBalazs et alJ 120031 ) . 
Some ha ve indicated the r e are more than tw o 
cluste rs dMukheriee et~aT| Il998t Irlorvathl Il998l ). 
while lHakkila et al] f|2000h maintain that the re- 
ported third cluster is simply the result of a se- 
lection bias. Furthermore, significant overlap is 
present in the duration comparison of short and 
long bursts, complicating a clear distinction be- 
tween two classes as has been discussed by the au- 
thors cited above. In any case, there are difficulties 



with all aforementioned classification measures in 
that the igo duration is somewhat subjective in the 
necessary selections of background regions, while 
a hardness ratio based upon counts is strongly de- 
tector dependent. 

The observation that some classically short 
bursts are extended in duration, when observed in 
an energy band (BATSE 20-50 keV) different from 
that of the n atural BATSE 50-300 keV band, was 
report ed by lLazzati. Ramirez-Ruiz, fc Ghisellini 
"200lh. The dedicated GRB mission, Swift 



(jBarthelmv et al ]l2005l) . has introduced additional 
issues, while reinforcing the e xtended emis s ion o f 

(2000) 



Norris et al 



short GRBs at lower energies, 
have introduced the time lag between broad en- 
ergy channels ('spectral lag') as a classifier: 'short' 
GRBs have approximately zero lag, while the 
'long' events have a lag that is significantly dif- 
ferent from zero. Indeed, having a near-zero lag 
is the basis for claiming that some bursts ob- 
served by the Swift BAT that are significantly 
longer than 2 seconds belong in the 'shor t' class 
(jNorris fc Bonnel]||2006l ; IZhang et alj|2006j) . 

Fol lowing on the analysis of Nakar fc Piranl 
(120051) on the Amati r elation (|Amati et al.ll2002j ). 
Band fc Preecel (|2005h investigated the implica- 



tions of combinations of several observable GRB 
parameters, derived from an extensive data set of 
GRBs observed by BATSE. The BATSE data set 
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used was a partial spectral catalog of the peak flux 



and fluence spectral parameters (Mall ozzi et al 



19951 ) . comp lete up to the end of the BATSE 



4B Catalog (jPaciesas et a.1.1 Il999h . The BATSE 



5B Spectral Catalog has now bee n completed 
([Goldstein. Preece fc Mallozzil l2010t) . which in- 
cludes all BATSE bursts with sufficient counts 
in the spectral data that they could be analyzed. 
Based on this comprehensive data set, we present 
a new GRB classification measure that is as di- 
agnostic as the tgo duration for classification, but 
does not rely on the subject ive choices require d 
for the durations calculation ( Koshut et al. X996). 



2. Motivation 

BATSE data has bee n used to study the -E pea k 
distributions of GRBs (jKaneko et al.l 120061 ) . and 
the time-integrated -E pe ak distribution for all 
bursts shows no evidence of discrimination be- 
tween short and long bursts. Fluen ce hardness 
distributions (Ko uveliotou et al. Ill993l) show some 
evidence of bimodality (jBalazs et al.l 120031) , but 
there is only mode rate significance with much 
overlap (|Nakarll2007l ) . Using the BATSE 5B Spec- 
tral Catalog, we support previous observations on 
the distributions of -E pea k and fluence. In addi- 
tion, we split the distributions into long and short 
duration GRB distributions, following the tgo clas- 
sification of two seconds. The Speak distribution 
for short bursts is completely overlapped by that 
of long bursts. Although the S pca k values for long 
bursts are centered around 150 keV, short burst 
Epcak values are shifted to higher energies around 
300 ke V. This is consisten t wit h previous find- 
i ngs o f IPaciesas et al. ( 2003 ) and Ghirlanda et al 



(2009). Approximately 65% of the fluence distri- 
bution for short bursts overlaps that of the fluence 
distribution for long bursts, with the position of 
the peak of the short GRB fluence distribution 
being an order of magnitude less than that for 

lo ng GRBs. 

Llovd. Petrosian. fc Mallozzil (2000) have shown 



there to be a significant correlation between S poa k 
and the total fluence in gamma-rays, and it is 
desirable to investigate this correlation for both 
long and short GRBs. For this reason, we investi- 
gate the Speak - fluence distribution for all BATSE 
bursts with good spectral fits and devise a discrim- 
inator between long and short bursts based on the 



difference in correlation between Speak and fluence 
for long and short bursts. A choice formulation 
for a d iscriminator is bas e d on the hardness of a 
burst. iKouveliotou et all (1993) used a hardness 
ratio based on the ratio of calculated fluence in 
different energy bands to compare to duration es- 
timates. Instead, we propose to use a hardness 
measure represented by S pe ak/fluence. This so- 
called energy ratio is in units of area and should 
prove to be a good discriminator between long 
and short bursts if there is a strong correlation 
between Speak and fluence. 

3. Observations 

From the 5B spectral catalog (Goldstein, Preece 
& Mallozzi, in prep.), we extract bursts with a 
good model fit as determined by a 3-sigma con- 
fidence limit, with the time-integrated S pea k and 
fluence errors for each burst required to be no more 
than 40% of their respective fitted values. A total 
of 1121 long bursts and 168 short bursts, classified 
according to the classica l tg cut of two seconds 



( Kouveliotou et al. 19931 ). satisfied these criteria. 



We then calculated the S pca k/fluence energy ra- 
tio for each of these bursts and plot a histogram 
of these values, as in Fig. 1. By using a stan- 
dard nonlinear least-squares fitting algorithm, we 
fit a single lognormal function to the distribution 
with the resulting chi-square goodness-of-fit statis- 
tic 111 for 32 degrees of freedom. We then fit two 
lognormal functions to the distribution with the 
resulting chi-square statistic of 32 for 29 degrees 
of freedom. Since the two models are nested, we 
use Pearson's chi-square test to show that the large 
change in chi-square per degree of freedom results 
in a chance probability of 5 x 10~ 17 and that the 
two lognormals are statistically preferred with a 
high degree of significance . 

From this bimodal distribution an obvious dis- 
tinction between long and short bursts emerges. 
In Fig. 2 we plot two histograms corresponding to 
long and short bursts as identified by their respec- 
tive tgo estimations to show that the bimodal dis- 
tribution of the energy ratio is correlated to that of 
the duration distribution. A K-S test comparing 
the long and short burst distributions in Fig. 2 to 
the best fit lognormal functions in Fig. 1, however, 
finds the correlations to be statistically marginal 
with 2% and 0.6% respective probabilities that 
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each distribution in Fig. 2 is drawn from their 
respective best fit lognormal functions in Fig. 1. 
The energy ratio distribution for short bursts is 
narrower compared to that of long bursts and is 
shifted to higher energies, resulting from the fact 
that short bursts are generally harder than long 
bursts. It appears the energy ratio values are a 
good discriminator between the classical definition 
of long and short bursts by merging two known 
discriminators into one quantity, and their rela- 
tive overlap can be well estimated. Only 4% of 
long GRBs overlap the 1-sigma core of the short 
burst distribution, and 2% of short bursts overlap 
the the 1-sigma core of the long burst distribu- 
tion. Similarly, the overlap for the 2-sigma cores 
is 11% and 23%, respectively. Comparatively, for 
our sample, the classical t 90 overlap of long bursts 
onto the 1-sigma (2-sigma) core of short bursts is 
4% (23%), and and there is no overlap of short 
bursts onto either the 1-sigma or 2-sigma core of 
long bursts. The central value for the long burst 
energy ratio distribution is ^0.06 while the central 
value for the short burst distribution is ~1.5. 



relation ( Amati et alj |2002) 



Band & Preecc (2005) showed that the Amati 
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lower limits. Using the bursts described above, 
the i? P eak values and fluences can be plotted in 
this plane. Fig. 3 shows the distribution of long 
bursts in the upper plot and the distribution of 
short bursts in the bottom plot. The lines denote 
the lower limits of the Amati and Ghirlanda rela- 
tions in this plane. Note that the beaming fraction 
for the Ghirlanda relation, /b = 1.0, which is re- 
lated to the jet opening angle, 9, by fs = 1 — cos 8. 
This represents the energy radiated by the burst 
equalling £% D . It can easily be seen there are two 
separate distributions in the Speak _ E- lso plane. 
Long bursts appear to be clumped between the 
Amati and Ghirlanda limits, while the short bursts 
appear to distribute along the Ghirlanda lower 
limit. 

4. Conclusions 

The energy ratio shows a clear distinction be- 
tween two different types of GRBs. The fluence 
encodes the duration of the burst without deriving 
a subjective tgo estimate, and -Epeak/fluence physi- 
cally represents a ratio of the energy at which most 
of the gamma-rays are emitted to the total energy 
emitted in gamma-rays. This quantity effectively 
serves as a spectral hardness ratio and shows an in- 
creased hardness for sh ort bursts compared to lon g 
bursts, consistent with iKouveliotou et all (|l993h . 



2004), 



and the Ghirlanda relation (IGhirlanda. Ghisellini fc 

(2) 



Lazzat 



n he distribution of short bursts is narrower, its 1- 



^peak-^G I 1055-^ 



could be converted into a similar energy ratio 
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<xF(z), 



(3) 



where the C\ in the previous equations are the re- 
spective normalization coefficients and /b is the 
beaming fraction relevant for the Ghirlanda rela- 
tion. Here, S 7 is the fluence in gamma-rays, and 
r\i are the best fit power law indices for the re- 
spective models. These energy ratios can be rep- 
resented as functions of redshift, F(z), and the 
upper limit of the ratios could be determined for 
any redshift. The energy ratio upper limit of the 
Amati relation, as well as the energy ratio upper 
limit of the Ghirlanda relation, can be projected 
into the S pea k ~~ fluence plane where they become 



sigma width covering less than one order of mag- 
nitude in energy, while the long bursts 1-sigma 
width covers slightly more than an order of magni- 
tude. This bimodal distribution heavily supports 
the ori ginal distinction betwee n long and short 
bursts (|Kouveliotou et al.lll993f ) and suggests fur- 
ther investigation is desirable. The correlation be- 
tween Speak and the total fluence in gamma-rays 
is of particular interest, as the energy ratio re- 
moves the cosmological dependence for the ener- 
gies involved since its value is merely proportional 
to the square of the luminosity distance. Because 
of this reason, the energy ratio could be consid- 
ered physically superior to the duration classifier 
for GRBs. The low degree of correlation between 
the tgo distribution and the energy ratio distribu- 
tion may be attributed to the fact that the former 
is a an observed quantity in the observer's frame, 
while the latter is a quantity that contains spectral 
information from the rest frame of the GRB. In ad- 
dition, when comparing the difference in overlap 
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between the energy ratio distributions and the tgo 
distributions, the overlap for the tgo distributions 
is marginally less pronounced, but gives little in- 
sight to the physical processes in the rest frame 
of the GRB. Note that there may be a truncation 
effect with the energy ratio associated with short 
bursts, due to the inability of BATSE to trigger 
on very low fluence events and bursts with -E pC ak 
outside the BATSE energy range. It is expected 
that Fermi/GBM will assist in discovering GRBs 
with -Epoak values greater than 1 MeV, alleviating 
possible truncation effects due to the energy cutoff 
at the high end of the detector energy band. 

While a majority of BATSE bursts (-87 per- 
cent) fail the lower limit test for the Amati rela- 
tion in the -Speak ~~ fluence plane, very few BATSE 
bursts violate the lower limit for the Ghirlanda 
relation with /b = 1.0. Especially intriguing is 
the fact that short GRBs fall close to the lower 
limit where the energy radiated is isotropic. De- 
creasing the beaming fraction shifts the lower limit 
to higher fluences, causing an increasing number 
of bursts to violate the Ghirlanda relation. In 
addition, only a few jet breaks have been dis - 
cove red for short GRBs ( Soderberg et al. 2006), 
and Watson et al. ( 20061 ) caution those discover- 
ies may be misleading due to the flaring activities 
of short GRB decay. Thus, the Ghirlanda rela- 
tion, as well as the lack of reliable observed jet 
breaks for short bursts, suggests that short bursts 
release energy isotropically (or near isotropically), 
while longer bursts tend to have a much smaller 
beaming fraction and consequently have small 
opening angles. Clearly the Ghirlanda relation, 
if accurate, requires short bursts to release en- 
ergy isotropically, as opp osed to beamed radiation 



2001; Nakar 



relea se in long bursts (jFrail et al 
2007). In addition, if short GRBs are isotropic 



emitters, then virtually all short GRBs should be 
detectable within a given volume of space. Long 
GRBs, in general, have a small measured open- 
ing angle, a nd therefore only a small fraction are 
detec table ([van Paradiis. Kouveliotou. fe Wiiers 
2000). Since short GRBs are detected far less fre- 
quently than long GRBs (jPaciesas et al]ll999j ) this 
is indicative of a relatively rare and independent 
cause for short GRBs. 
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Fig. 1. — Histogram of the energy ratio distribu- 
tion in the 20-2000 keV range for 1289 GRBs. The 
dotted line shows the best fit lognormal to the dis- 
tribution, and the dashed lines show the fit of two 
lognormal functions. The chi-square goodness-of- 
fit statistic for one lognormal is 111 for 32 degrees 
of freedom, while for two lognormals is 32 for 29 
degrees of freedom. Therefore, two lognormal dis- 
tributions are statistically perferred, resulting in a 
bimodal distribution for the energy ratio. 



6 



100 




0.001 0.010 _ 0.100 1.000 10.000 

Epeak/Fluence 

Fig. 2. — Histograms of the energy ratio distri- 
butions in the 20-2000 keV range for 1121 long 
bursts (white) and 168 short bursts (gray). There 
are clearly two distinct distributions, with long 
bursts centered around 0.6 and short bursts cen- 
tered around 1.5. The solid curves are the best fit 
lognormal functions, and the dashed lines are the 
1-sigma standard deviation of the distributions. 
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Fig. 3. — Scatter plots of 1121 long bursts (top) 
and 168 short bursts (bottom) in the i? poa k 
fluence plane. Also plotted are the lower lim- 
its of the Amati relation and Ghirlanda relation 
(fB=1.0). The dashed lines represent the 1-sigma 
errors about the lower limits. There is apparently 
a clear lower limit violation by most BATSE bursts 
for the Amati relation. Short bursts appear to 
cluster around the lower limit of the Ghirlanda 
relation. 



